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Abstract: Dark matter direct detection searches for signals coming from dark matter scat- 
tering against nuclei at a very low recoil energy scale ~ 10 keV. In this paper, a simple 
non-relativistic effective theory is constructed to describe interactions between dark matter 
and nuclei without referring to any underlying high energy models. It contains the mini- 
mal set of operators that will be tested by direct detection. The effective theory approach 
highlights the set of distinguishable recoil spectra that could arise from different theoretical 
models. If dark matter is discovered in the near future in direct detection experiments, a 
measurement of the shape of the recoil spectrum will provide valuable information on the 
underlying dynamics. We bound the coefficients of the operators in our non-relativistic effec- 
tive theory by the null results of current dark matter direct detection experiments. We also 
discuss the mapping between the non-relativistic effective theory and field theory models or 
operators, including aspects of the matching of quark and gluon operators to nuclear form 
factors. 
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1. Introduction 

The existence of dark matter (DM) in the Universe is well established by astrophysical and 
cosmological observations. However, these observations only detect DM through the univer- 
sal interaction of gravity, leaving the nature of DM and its interactions almost completely 
unknown. Attempts to understand DM more directly by looking for signals of its interaction 
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with the Standard Model (SM) are taking place on a variety of experimental and observational 
frontiers. Among them, direct detection searches for signals from DM scattering off nuclei 
in underground detectors. In these experiments, signals have traditionally been assumed to 
come from the simplest form of elastic scattering, a contact interaction independent of mo- 
mentum exchange. This assumption leads to the constraints from the present null results of 
most direct detection experiments, usually presented as bounds on the spin independent (SI) 
and spin dependent (SD) cross sections per nucleon as a function of DM mass. 

The simple assumption of a DM contact interaction with nucleons overlooks the possible 
sensitivity of direct detection to more general DM scenarios. This has already been pointed 
out in the context of inelastic DM [1, 2] and form factor DM [3, 4]. In general, the form of the 
DM-nucleus scattering amplitude depends on properties of the dark matter particle as well 
as the mediator of the interaction. As we will show in this paper, other interactions could 
also generate sizable direct detection signals. Considering the large number of existing DM 
models or effective field theory operators for DM with different spins, it may appear difficult 
to exhaustively categorize all possible interactions that could be tested by direct detection. 
However, scattering in direct detection is non-relativistic (NR) as the incoming DM velocity 
is v/c ~ 10~ 3 . Thus the incident DM kinetic energy and recoil energy are around 10 keV. 
At this low energy scale probed by direct detection, many apparently different microscopic 
models or field theory operators lead to the same simple NR effective theory. In other words, 
the DM-nucleus scattering in direct detection can be essentially described by an NR effective 
potential with small expansion parameters: the DM velocity v ~ and |<f|/A. \q\ ~ O(10 
- 100 MeV) is the momentum exchange and A is some large scale involved, such as the DM 
mass m x , the nucleus mass m^, or a heavy mediator mass. 

This NR effective theory description serves as a systematic approach to parameterize the 
direct detection signals for various underlying microscopic DM theories. It captures the most 
important differences in the measured recoil spectrum for different classes of DM models. 
Thus it is the simplest theory with the minimal set of operators that will be tested by the 
direct detection. Current direct detection results bound the coefficients of operators in the 
NR effective potential. Constraints on specific DM models or field theory operators could be 
obtained by mapping them onto the NR effective theory. 

This paper is organized as follows: in section 2, we present the NR effective theory for 
direct detection and complete the various NR operators into microscopic models. In section 3, 
we apply current direct detection data to constrain the Wilson coefficients of our NR theory; 
we also discuss constraints from other experiments on different force carriers between DM and 
nucleus. In section 4, we discuss the procedure of mapping from relativistic quark or gluon 
operators to NR nucleus operators. We conclude in section 5. Finally, we include formulae 
for the nuclear recoil rate, power counting rules and the mapping between NR theory and 
field theory operators in appendices. 
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2. NR effective theory for direct detection 
2.1 General considerations 

The DM velocity in the galactic halo is typically of the order v ~ 10 -3 . For DM and 
target nuclei with mass around the weak scale, this implies incident DM kinetic energy and 
recoil energy around 10 keV, much smaller than typical nuclear binding energy (1 - 10) MeV 
per nucleon. This justifies the treatment of DM scattering against the whole nucleus using 
NR quantum mechanics. In our framework, the DM-nucleus interaction is described by an 
effective potential, V e g. It is a rotationally invariant scalar formed out of four 3-vectors: the 
relative velocity v, DM position r, DM spin s x and nuclear spin sjy. The case in which sn 
does not appear is known as spin-independent (SI) scattering, as opposed to spin-dependent 
(SD) scattering. To the first order of the Born approximation, the amplitude in the NR limit 
is the Fourier transform of the effective potential in position space: 

M(q, v) = ~J d 3 re^ r X s (r, v), (2.1) 

where q is the transferred momentum. In general, the f dependence will involve the structure 
of the nucleus, due to its finite size. This effect requires a nuclear form factor depending on 
the transferred momentum q in the calculation of the differential recoil rate. In tabulating 
the operators that can appear in V e s, we will factor out AF{q), the number of nucleons times 
the nuclear form factor, in the SI case, and y/S{q)/S(0), the spin form factor, in the SD case. 
(So, for instance, if the DM interacts with the charge of the nucleus, we factor out AF{q) and 
the coupling will be rescaled by a factor of Z/A.) We will discuss more details of the role of 
nuclear physics and the form factors that can appear in Sec. 4, and give explicit examples of 
the matching procedure in detail in App. C.1.3 to clarify our conventions. 

The mass scales that enter into the potential are the DM mass m x , the nucleus mass mjv, 
and the mediator mass m^. In addition, there could be other scales A' present. For instance, 
the DM-mediator interaction could arise at the nonrenormalizable level, encoded by a high 
dimensional operator suppressed by powers of A'. This could happen, for example, in models 
of DM with zero electric charge but higher-order electromagnetic form factors. Notice that 
direct detection experiments constrain one combination of the mass scales and the coupling 
constants. Thus, in the following parameterization of NR effective theory, we choose to absorb 
all of the scale dependence into the coefficients of the terms in the effective potential. 

We will consider two qualitatively different cases: a contact interaction when |g| <C 
and the mediator is integrated out, and a long range interaction when |g| 3> m^. The 
leading term in the effective potential (after factoring out the nuclear form factor) is then 
proportional to 5 3 (f) and 1/r, respectively. For a contact-type interaction, the NR expansion 
of the potential is a derivative expansion, while for a long-range interaction, it is a multipole 
expansion. It is often assumed that, due to the smallness of the expansion parameters, direct 
detection experiments will only be sensitive to a momentum-independent potential. However, 
as will be shown in Sec. 2.3, it is possible that the leading contributions from microscopic 
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models are already suppressed by a single |g| or \v\. High-order terms in \q\ could still be 
detectable if they are leading, especially for the case of long range interactions. However, given 
the parameter space probed by experiments, the potentials relevant for direct detection form 
a finite set. In the studies below, for simplicity, we will focus on potential terms suppressed 
by at most a single | q\ . More details of power counting rules for the two cases may be found 
in Appendix B. 

2.2 NR effective potential 

We begin by writing down the effective potentials for SI scattering and V^ D for SD 
scattering, 

v cS = vM + v c f 

= M 3 (^ - h 2 s x ■ W 3 (f) 

+ h-— + h-^— (2 2) 

47rr 4-7rr 3 ' 

KfP = h[s x ■ s N 5 3 (r) - h' 2 s N • W 3 (r) 
1 4nr 2 4nr 3 

where the (dimensionful) Wilson coefficients h (for heavy mediators) and I (for light medi- 
ators) are determined by matching the field theory operators from the underlying theory to 
the NR operators. They are proportional to the couplings of DM to the mediator, g x , as well 
as the couplings of the nucleus to the mediator (which in general can involve nuclear physics 
quantities like those traditionally denoted in the DM literature). One should understand 
the terms in the effective potential to carry spin indices, which we have suppressed; for in- 
stance, h\5 3 {r) multiplies 5 rr t5 ss i, where r,r' and s,s' are the spins of the nucleus and DM 
before and after scattering. Similarly, s x is to be understood as an appropriate representation 
of spin, s x - ss /] e.g. for fermionic DM, it will be given by the Pauli matrices in a basis with 
given spins, ^^ s '^a^ s \ Numerical factors in the translation of V c s to a spin-averaged cross 
section da/dEn will depend on the spin structure. We will adopt the convention that all cross 
sections and limit plots are for the case that DM is a Dirac fermion; the results can be easily 
rescaled to other cases. Detailed examples of matching are given in section C.1.3. Note that in 
general, the couplings h[, . . . will be isotope-dependent. DM models that give rise to the above 
NR operators have already appeared in the literature (Recent examples include momentum 
dependent DM [3, 4] and DM with electromagnetic form factor [5, 6, 7, 8, 9, 10, 11, 12, 13]). 
The NR effective theory will manifest the experimentally testable physics without referring 
to diverse high energy interpretations. 

The first operators in and V^ D are the most-studied cases, SI and SD contact inter- 
actions. The momentum suppressed operators are usually neglected. Yet in the absence of 
momentum-independent operators or with large enough coefficients to compensate the mo- 
mentum suppression, these terms could still be relevant for direct detection. A couple of 
comments are in order: 
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• Any specific DM model typically gives a subset of the operators presented in Eq. (2.2) 
and Eq. (2.3). The most stringent bound from direct detection is always only on the 
leading order operator. As we will show, each operator presented is the leading contri- 
bution from a natural UV completion. 

• We have assumed a derivative expansion in r, corresponding to powers of q in momentum 
space. Thus, for instance, a 1/r 3 term is omitted because it corresponds to a log|tf| 
interaction. Logarithmic dependence on \q\ signals the exchange of a continuum of 
modes, and could arise if both DM and the SM couple to a new massless or conformal 
sector (which could have a mass scale far below the other scales in the problem). Such 
sectors have received attention in the guises of RS2 [14] and unparticles [15], but we are 
not aware of models of such a sector mediating interactions between DM and the SM. 
A log | q\ dependence would have fairly mild effects on the shape of the recoil energy 
distribution, which is subject to uncertainties including nuclear physics and the DM 
velocity distribution, so we expect that it would be extremely challenging to deduce 
evidence for a new conformal sector solely from DM direct detection. 

• We consider only static potentials. There are also operators depending on the relative 
velocity v of the DM and nucleus. For instance, the operators suppressed by a single v 



It turns out that for elastic scattering, numerically the shapes of the recoil spectra 
they generate are almost identical to those of the unsuppressed operators. (Their ef- 
fect is most important on the tail of the Er distribution, which is relatively unimpor- 
tant.) Thus we neglect them for the rest of the studies. Experimentally, this higher- 
dimensional space of potentials leads to a "degeneracy" in the inverse problem of map- 
ping measured spectra to NR operators. We discuss this in somewhat more detail in 
Appendix C.1.3. 

• In eqs. (2.2, 2.3), we only keep operators suppressed by a single \q\. Our potential is 
the minimal extension of the well-studied momentum independent contact interaction. 
If the leading operator from a field theory is suppressed by more powers of | q\ , it could 
still contribute an observable rate to direct detection, depending on the mediator mass. 
We will give a more exhaustive list in Appendices C.l, C.2, C.3 but focus on eqs. (2.2, 
2.3) in the numerical studies. 

As emphasized in the introduction, the effective theory parameterization highlights the 
possibility of having qualitatively different recoil energy spectra, as shown in Table 1. 

2.3 Connection with microscopic models 

The matching between microscopic models and effective potentials goes through three straight- 
forward steps. Starting with a particular DM model, we first write down the relativistic field 
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SI NR operators 


SD NR operators 


Er 


5 3 (r) 

s x • W 3 (f) 
1 

4nr 
4irr 3 


s x ■ s N 5 3 (f) 
s N ■ W 3 (r) 

s x -s N 
4irr 

4 7rr 3 


1 

Er 
Er 2 
Er 1 



Table 1: Recoil energy dependence of effective cross section per nucleon for the NR operators. (Nu- 
clear form factors are factored out.) 

theory operators relevant for the scattering process, G(q 2 ,m ( p)J x J q ^ g y J x and J q ( g ) are ap- 
propriate DM and quark (gluon) operators, respectively. G(q 2 ,m ( f ) ) comes from the exchange 
of the mediator. Second, we convert this to an operator involving nucleons by taking the 
nuclear matrix element (N(p + q)\ J q ( g \(q)\N (p)} . Then we take the Fourier transform of the 
scattering amplitude in the NR limit to get the effective potential, factoring out the nuclear 
form factor in the definition of V e &. We present the result of the final step matching in the 
cases of scalar, fermion and vector DM in Appendix C.l, C.2, C.3. Details of the nuclear 
matrix elements are discussed in section 4. In this section, we discuss simple high energy 
models of fermion DM for the set of NR operators in Eq. (2.2, 2.3). 

Most of the existing DM models yield simple contact interactions as the leading operator. 
It is possible, however, the coefficients of momentum (velocity) suppressed operators dominate 
over that of the simplest contact interaction. The coefficient enhancement could be due to 
large couplings of a mediator to the dark sector or small mediator mass as a ~ m~ 4 . In 
models where the leading operator is SD while the SI operators are momentum-suppressed, 
the SI scattering could still be detectable as SI searches probe much weaker processes than 
SD searches. For the SI operators, we have 

• 5 s (?) 

This is the most studied case. For instance, Higgs exchange between fermion DM and 
the nucleus would lead to an operator XXQQ- With Higgs mass around 100 GeV and 
the Higgs-nucleon coupling about the strange-quark Yukawa coupling 10~ 3 , this would 
lead to a plausible scattering cross section ~ 10 -44 cm 2 . This could also arise from 
gauge boson exchange, e.g., a Z boson exchange. The DM-Z coupling has to be of 
order 10~ 3 for a 10~ 44 cm 2 cross section. This small coupling could come from some 
high dimensional operator. For instance, assume that DM and the SM higgs are both 
charged under a new U(l) which is broken at scale A ~10 TeV. Integrating out the Z' 
leads to a dimension six operator xi^xh D^h which induces an effective DM-Z coupling 
of order (u w /A) 2 ~ (10~ 4 - 10~ 3 ). 

. s x -W 3 (r) 

This operator could arise from d'x^^l^xE 1 ^, an electric dipole coupling of DM to a 
new GeV gauge boson which kinetically mixes with the photon eF' F^ u [16]. The NR 
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operator coefficient hi is related to the dark electric dipole moment d! as hi ~ d'e/m\. 
The direct detection bound hi < 10~ 7 GeV~ 3 becomes die < 10~ 21 (e cm). Notice that 
this interaction leads to a potential proportional to the charge Z of the nucleus instead 
of the atomic number A. 

This operator could also be generated by ixi^XQQ horn heavy scalar/pseudoscalar ex- 
change. In this case, hi ~ gj (m x m'^ ) ), where the dimensionless coefficient g is propor- 
tional to the CP-odd coupling from either the DM sector or the visible sector. For 
DM with mass around the weak scale m x ~ 0(100 GeV), g ~ O(0. 1) for a 100 GeV 
mediator and g ~ O(10 -5 ) for a 1 GeV mediator. 

This operator could also appear in a special linear combination with the velocity- 
suppressed operator s x ■ v, giving a shape that is similar to a linear combination of the 
shapes from couplings h\ and hi. The example is the DM "anapole" moment scattering 
off the nuclear electric current, which will be discussed in detail in Appendix C.1.3. 

• 1/r 

The Coulomb potential could arise through exchange of some new light bosons with 
< \q\. The current direct detection bounds require the coupling to be tiny, ^ 
10 -11 (see Sec. 3.2). Although such a coupling may seem unreasonably small, it is 
not difficult to build dark photon models that satisfy the bound, e.g. the new light 
gauge boson kinetically mixed with the photon only through S'-parameter-like higher 
dimension operators, as discussed in [17]. 

• s x • r/ r 3 

This is the DM-dipole/nucleus-monopole coupling. They could arise from models similar 
to those that generate the contact interaction s • V5 3 (r), but with the mediators lighter 
than the momentum transfer. For instance, DM electric dipole coupling to the ordinary 
photon would give this potential with coefficient li ~ d. The direct detection bound 
li < 10 -9 GeV -1 gives the DM electric dipole moment d < 10 -23 (e cm). For models 
with light scalar or pseudoscalar exchange, li ~ g/n^x lading to the CP odd coupling 
g < 10 -7 for weak scale DM. 

Now we would like to estimate roughly how large the CP violating coupling could be 
if the CP violation is confined to the visible sector. CP-violating phases in the visible 
sector could be present in various extensions of the SM. They generically have to be 
small to avoid generating electric dipole moments for the neutron, electron, and atoms 
in conflict with observed data. For instance, in the minimal supersymmetric SM, a 
combination of complex phases of the gaugino-mass parameters, the A parameters, and 
fi must be less than the order of 10 -2 — 10~ 3 (for a supersymmetry-breaking scale 
of 100 GeV). If a CP violating coupling like vr'qg with ir' a pseudoscalar is generated 
at one loop, the coupling may be further suppressed by the loop factor. Taking into 
account other possible small couplings, one could get values in the range 10 -5 — 10 . 
Intriguingly, this is around the bound set by direct detection. 
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For the SD operators, the leading contact interaction s x ■ sjv5 3 (r) again arises from most 
existing models, e.g, from Z exchange or squark exchange in the neutralino DM scenario. 
If the mediator gauge boson is light, it may lead to the long range SD potential s x ■ sjy/r. 
Consider the case that the mediator is a pseudoscalar with CP- violating coupling to the dark 
sector tt'xx due to some hidden CP violating phases. If the mediator is heavy, it will lead 
to Sjv • V(5 3 (r), while if it is light, it generates the nucleus-dipole/DM-monopole potential 
sn • r/r 3 . 

3. Constraints and sensitivities 

3.1 Constraints from current direct detection experiments 

In this section we constrain the coefficients appearing in Eq. (2.2, 2.3) using the most sensitive 
direct detection experiments. We consider two classes of experiments based on different 
techniques. One class of experiments is based on measurement of the energy DM deposits in 
a detector by scattering off target nuclei. Most direct detection experiments, such as CDMS- 
II [18, 19], XENON10 [20, 21, 22], and XENON100 [23], belong to this class. Another class 
of experiments, especially useful for the SD constraints, is based on the superheated droplet 
(bubble chamber) technique to search for DM recoiling on 19 F nuclei in a compound target. 
This process is very sensitive to SD interactions. If the energy deposited by DM exceeds a 
minimal energy barrier determined by the experimental thermodynamical conditions, it will 
trigger the nucleation of a bubble of the gas phase in the superheated liquid. The experiments 
search for bubble formation events not accounted for by various backgrounds. They constrain 
the total rate of DM-nucleus scattering but cannot give the differential rate as the first class 
of experiments do. COUPP [24], SIMPLE [25] and PICASSO [26] are the leading experiments 
in this class. 

For the SI constraints, we use data from CDMS-II [18, 19], XENON10 [20, 21], and the 
currently available results of XENON100 [23]. For the SD case, we use data from XENON10 
and XENON100 for the DM-neutron coupling, as half of naturally occurring xenon is in the 
form of isotopes with unpaired neutrons [22]. We use the preliminary results of COUPP 
[24] to constrain the SD DM-proton coupling. For CDMS-II, XENON10, and XENON100, 
we will use the maximum gap method [27] to set limits while for COUPP, we will take the 
90% CL Poisson limit from their three DM candidate events, allowing ~ 6 events. In our 
calculations we have used form factors from Refs. [28, 29]. 

We use the parametrization of Eq. (2.2, 2.3) for numerical studies. The formulae for recoil 
rates are presented in Appendix A. We set one coefficient to be nonzero at a time and plot 
the spectrum in Fig. 1. From the left panel of Fig. 1, one can easily see that the cross section 
scaling as Er peaks toward high energy and is broader than the other distributions. The 
other rates peak at the threshold and have different slopes determined by their Er power 
dependence. In Fig. 1, we also show a spectrum with two comparable contributions from 
different operators. The spectrum resembles that of the "semielastic DM" scenario [30], as 
it peaks at the threshold but also rises at higher energies. The only difference is that in 
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"semielastic DM" the spectrum comes from a combination of elastic and inelastic scattering, 
while here we only consider elastic scattering. One could play further and consider both 
elastic and inelastic scattering and various combinations of operators. If, in the near future, 
direct detection confirms DM signals and collects enough data, one in principle could fit 
the recoil spectrum to our NR theory parameterization to obtain information about the DM 
interaction. 
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Figure 1: Left: Recoil energy spectrum for 150 GeV DM scattering on germanium. The curves 
correspond to spectra with one of hi,]i2,h,h being nonzero in red, blue dashed, purple dotted, and 
green dot-dashed. Right: Recoil energy of 150 GeV DM on germanium with nonzero /12 and l\. All 
the curves are normalized to have the same number of total events above 1 keV. 



The limits on each coefficient from the experiments we consider are displayed in Fig. 2 
and Fig. 3. Notice that in cases with a light mediator, due to the enhanced rate at small recoil 
energy, the XENON10 reanalysis from Ref. [21] becomes the strongest constraint, as it is the 
result with the lowest energy threshold. So far the importance of energy calibration around 
the threshold has been mostly emphasized for ruling in or out the light DM scenario [31]. 
For Xenon-based experiments, the recoil energy threshold is determined by the scintillation 
efficiency £> e ff: which is still controversial at the moment [23, 31]. But as is clear from the 
discussions above, the low energy region and the precise measurement of C e ff are also crucial 
for constraining heavy DM with a light mediator. From Fig. 2, the strongest constraints on 
the coefficients of SI operators in Eq. (2.2) are 
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(3.1) 



We have written the couplings as a ratio A, where the dimensionless part g ~ g x fp(n) is 
a product of the strength of the DM coupling to the mediator and the nucleus coupling 
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Figure 2: SI exclusion curves for Wilson coefficients with CDMS black dashed, XenonlO green solid, 
XenonlOO purple dotted. 



to the mediator, and we have chosen the scale involved, A, to be 100 GeV for illustrative 
purposes. (Parametrically, A could be m x ,mN,mj), or some new scale in the dark sector). 
From this decomposition, one can draw clues about the relevant scales and couplings of a 
model that could give rise to the NR operators. For instance, as already discussed in Sec. 2.3, 
the simplest contact interaction could come from Higgs exchange between DM and nucleus 
with Higgs mass around 100 GeV and its coupling to nucleus O(10 -3 — 10~ 4 ). We also include 
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in Fig. 2 an exclusion plot in (h,2, h) plane for 50 GeV DM. 

We plot constraints on the SD operators in Fig. 3. We have assumed purely neutron or 
purely proton couplings, for simplicity. In the case of the DM-neutron coupling, bounded 
by XENON10 and XENON100, the analysis is complicated by the fact that there are two 
important isotopes, 129 Xe and 131 Xe, with different spin interactions. We assume that these 
couplings scale as universal numbers h'i , , l'^ , l 2 times the isotope-dependent factor 

(S n ) /J- 



h[( A Xe) = h'} n) f^f) 



, (A = 129, 131) (3.2) 

and similarly for the other couplings. For 129 Xe, J = 1/2 and we take (S n ) = 0.359; for 
131 Xe, J = 3/2 and we take (S n ) = —0.236 [22]. We plot limits on h'j n \ etc., rather than on 
the isotope-dependent couplings bl x . The DM-neutron coupling strongest bounds are 

10~ 2 



h{ ' & 10" b GeV- 



2 



(100 GeV) 

h'^ & 10- 4 GeV~ 3 = 10-1 , 

(10 GeV) 3 

l[ {n) <, io- 8 

^^0-GeV--^^. (3.3) 

In the DM-proton case, no similar ambiguity arises, as we assume the scattering at COUPP is 
entirely on protons in the single isotope 19 F. Nonetheless, in order to have roughly comparable 
bounds in the two cases, we define h'^ by bl x = h'i {{S p } /J)|i 9p , and similarly for the 
other couplings. The DM-proton coupling constraints are a little bit weaker than the DM- 
neutron constraints. Notice that for the momentum-suppressed SD contact interaction to be 
detectable, it is necessary to have either a large coupling or the scales involved have to be 
small < O (10 GeV). 

3.2 Constraints on force carriers from other experiments 

In this section we briefly comment on constraints from other experiments on the mediators. 
Many constraints are model-dependent and mild compared to the direct detection constraints. 
A complete discussion is beyond the scope of this paper. 

First we consider constraints on heavy mediators from colliders. The interaction between 
the DM and quarks and gluons will also lead to direct DM production at hadron colliders such 
as the Tevatron. With initial state radiation, this may lead to mono-jet signals. Recently 
there were two studies of constraints on the contact interaction between DM and colored 
particles from monojet + missing Et searches at the Tevatron [32, 33]. They have shown 
that for a heavy mediator with > 100 GeV and SI interaction, collider searches have 
greater reach a ~ 10 -40 cm 2 for light DM with m x < 10 GeV, to which direct detection 
experiments have limited sensitivity. But for m x > 10 GeV, direct detection has a more 
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stringent bound. For SD interactions, the situation is more interesting. The Tevatron bound 
on cross section of contact interaction is (10 -39 — 10 -40 ) cm 2 for m x < 200 GeV, one or 
two orders of magnitude below that from direct detection. But once the mediator mass 
drops below 100 GeV, the production process is no longer from a contact interaction in the 
collider. The collider constraints weaken and the direct detection experiments have stronger 
or comparable bounds for m x > 10 GeV for both SI and SD interactions. 

For a light mediator with mass around or below a GeV, there could be model dependent 
constraints on its coupling to quarks from rare meson decays. Depending on the flavor 
structure of the couplings, rare Upsilon and kaon decays could be relevant. For example, 
T — > 70 constrains the couplings to the third generation of quarks and hence couplings to 
gluons (see Sec. 4.1 for the relation between couplings to heavy quarks and gluons). If the 
mediator couplings to the first two generations of quarks are not flavor diagonal, e.g. if l\ds(j) 
exists, there will be constraints from K + — > 7r + 0. Precise bounds depend on the final state 
of 4> decays. For instance, if 4> — > invisible, the experimental bound on the kaon branching 
ratio Br(K + — > 7r + + invisible) < 7.3 x 10 -11 gives l\ < 10~ 12 [34]. 1 More constraints from 
meson decays for vector mediators can be found in [35]. 

While mediator couplings to leptons are free parameters in general, if they are nonzero, 
they will be constrained by g — 2 measurement, beam-dump experiments and astrophysical 
constraints. First, the g— 2 experiment rules out large couplings to the leptons above 10~ 3 . For 
gauge boson mass below m„, da^ from axial vector coupling has an additional enhancement 
factor from the propagator, m^/m^, compared to 5a^ from other bosons and couplings [36]. 
In this case, the bound gets stronger with the coupling c < 10 _6 m^/MeV. The beam-dump 
experiments exclude couplings of order 10 -3 — 10™ 7 depending on the mediator mass. There 
are still small allowed islands in the parameter space around ~ (10 - 100) MeV and c ~ 
(10~ 5 — 10~ 7 ). Detailed analysis could be found in [37, 38]. For even smaller couplings, the 
constraints are from supernova cooling. If the light particle decays > 10 km from the point 
of production inside the supernova, it will contribute to the energy loss. To be compatible 
with the observed SN 1987 pulse duration, the new channel's contribution to the energy loss 
is estimated to be smaller than 10 53 erg/s. For vector mediator with mass below 100 MeV, 
this requirement rules out couplings around 10 -6 — 10 -8 [37]. For axion-like mediator with 

~ (1 - 200) MeV, the excluded coupling is of order 10 -9 [39]. 

4. Nuclear physics and form factors 

There is a subtlety involved in describing DM direct detection in terms of an effective field 
theory. The momentum transfers involved are q ~ 100 MeV, so 1/q can be of the same order 
as nuclear radii R ~ A 1 / 3 fm (recall that 1 fm _1 ~ 200 MeV). This implies that the nuclei 
are not pointlike from the perspective of DM, and we should not, strictly speaking, describe 
DM-nucleus interactions in terms of contact interactions V(r) ~ <5 3 (r); the DM-nucleus cross 
section will involve a form factor \F(q 2 )\ 2 . 

J Pure pseudoscalar coupling does not mediator this process unless it, <f> mixes. 
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v esc = 500 km/sec, p = 0.3 GeV/cm 3 w csc = 500 km/sec, p = 0.3 GeV/cm 3 
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Figure 3: SD exclusion curves for Wilson coefficients with Xenon 10 green solid, XenonlOO purple 
dotted (DM-neutron), COUPP black dashed (DM-proton). 
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Other aspects of nuclear physics can also influence the treatment of NR effective field 
theories. We have three small parameters, vdMiQ/M x , and q/Mjy, all of which are of order 
10 -3 . However, nuclear binding energies are of order 10 MeV per nucleon, so inside the 
nucleus, individual nucleons have velocities v n ~ 10 _1 . One might worry that this, along 
with q/M n , is a large expansion parameter that could spoil the NR effective theory. We will 
review the treatment of nuclear form factors in DM scattering to argue that this does not 
occur. 

We would also, in the spirit of effective field theory, like to give an overview of the various 
operators that can be involved in nuclear scattering. A full review of those aspects of the 
physics of nucleons and nuclei that are relevant for DM direct detection is beyond the scope of 
this paper. Nonetheless, we will point out some lacunae in the DM literature. In some cases, 
we will offer pointers into the QCD and nuclear literature that can begin to fill these gaps. 
The early direct detection literature understandably focused on the operators q^^bq and qq 
that are most relevant for direct detection in the MSSM, linking their matrix elements to 
the nucleon spin and the pion- nucleon sigma term [40, 41, 42]. The nuclear physics involved 
in these operators has been reviewed many times [43, 44, 45]. More recent literature has 
broadened the set of operators considered, but often a discussion of nuclear matrix elements 
has been omitted. 

4.1 Matching to relativistic operators 

First, let us discuss the most basic step of the matching, from a renormalizable QFT to 
effective operators of the form J q ^J x . Suppose we integrate out a heavy mediator, say a 
scalar of mass M exchanged between DM and a quark. The leading term will be a contact 
interaction XXQQ f° r which we care about the form factor of a scalar current qq in the nucleus. 
What about subleading terms? The full tree- level amplitude is XX 2_ 1 A /2 gg 1 which has an 

2 

expansion in 4^ but in which every term will depend on the same form factor (N\qq\N). 
When integrating out the heavy scalar, one could also draw diagrams with radiated gluons, 
but this generates operators suppressed by both factors of a s (M) and further powers of M, 
and we can safely neglect them. (In certain cases one may wish to include corrections involving 

(2) 

derivatives, as in the twist-two coupling Xl^^X^q^v studied in [46]. In any case, one expects 
that in general operators with many derivatives are unimportant.) 

Next one can consider again the case of exchange of a particle between DM and a quark, 
but in the case that this particle is light, so that the amplitude is XX-^qq- Again the form 
factor is (N\qq\N). Unlike our previous argument, one cannot now argue that radiating 
additional gluons is suppressed by small a s or powers of the mediator mass M. However, one 
doesn't have to: in this case, we are not integrating out a particle and matching onto local 
operators at a high scale M, so we should never think of the quarks as perturbative. We 
simply have a process with amplitude xx^zQQ m t ne low-energy theory, and we are evaluating 
the matrix element (N\qq\N) in the nonperturbative QCD regime. Any gluons radiated from 
the quarks should be thought of as part of the nonperturbative "blob" connecting the local 
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operator qq to the nucleoli external states. 

Another possibility is that we have an s-channel particle, like a squark in a process 
XQ — > q* — > XQ- I n this case the amplitude has the form XQ s J~m 2 so the q and q are not at 
the same point and in principle we care about a form factor for a nonlocal operator, which 
might be thought of as involving contributions of infinitely many local operators. However, 
we are saved by the fact that s ~ (tun + m x ) 2 3> (q 2 , Aq CD ), so that for all practical purposes 
the operator we generate may again be treated as local (up to g 2 /s-suppressed terms). We 
can Fierz (xq)(xq) to obtain an expression in terms of local gauge-invariant operators like 
qq and qj p q, and then use the simple form factors {N\qTq\N) . Note that this conclusion is 
independent of the mediator mass M, provided the process is not resonant with s ~ M 2 . In 
particular, even though one cannot strictly "integrate out" a light s-channel particle, for the 
purpose of understanding nuclear form factors we can still treat this case as leading to a local 
operator qq. 

DM can interact not only with quarks but with glue. The simplest way to generate an 
operator like XXG^ lu G aiJ ' u is if DM exchanges a scalar (like the Higgs) with heavy quarks. 
Alternatively, DM could exchange a pseudoscalar with heavy quarks. In both cases we can 
integrate out the heavy quarks to produce effective couplings to gluons. The matching is [47]: 

m Q QQ -> -\~G%G a » V (4-1) 

m Q Q l5 Q -> i^G^G^. (4.2) 

It is natural to then ask what happens when integrating out heavy quarks exchanging other 
types of mediators. This possibility has received relatively little attention in the DM context, 
but most of the ingredients necessary to understand it exist in the QCD literature. The back- 
ground field method gives an efficient technique for computing the effective gluonic operators 
generated in such cases [48]. For instance, a new vector coupling to a heavy-quark current like 
Q^Q must match to some conserved current below the scale mg. One generally expects that 
conserved currents are dimension 3, but there is another possibility: they can be descendants 
of other operators, rather than primaries, so that matches on to an operator of the form 
du T )jLv with t^u an antisymmetric tensor. This must be the case in order to match to a gluonic 
operator after integrating out the heavy quark. The result is [49]: 

Q^Q -)• (167r2) f 45m 4 (?>d a TrG aT {G ar , G afl } - Ud a TrG^{G° T , G Ta }) (4.3) 
QWQ -> 48 jt 2 ^rraTr (G ap d a G T<7 + 2G Ta d a G ap ) . (4.4) 



m 



On the other hand, if we have a dipole-type interaction with Q^cr^Q, then integrating out 
the heavy quark generates couplings to dimension-six gluonic antisymmetric tensor operators. 
Again the matching is best calculated using the background field method [50, 51]: 

QT6<VQ -> — 3- d abc G a G hpa G c . (4.5) 
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These operators are suppressed by powers of the heavy-quark scale, so we expect that they 
would be overwhelmed by contributions from light-quark currents except in models where DM 
interacts primarily with the third generation through new forces. It would be interesting, but 
beyond the scope of this paper, to consider how the size of effects in direct detection and 
flavor physics might be related in such models. 

These arguments convince us that, despite the in- principle issue that DM-nucleus scat- 
tering can depend on arbitrary form factors {N(p + q)\(D(q)\N(p)), in practice only a handful 
of form factors (for operators made up of two quarks or two or three gluonic field strengths) 
will be relevant and contributions of others will be suppressed by powers like q 2 /M 2 or q 2 /m 2 N , 
sometimes with small perturbative couplings a s in front. 

4.2 Matching quark and gluon operators to nucleon operators 

Once we have an effective theory involving DM currents coupled to quark or gluon currents, 
the next step is to match onto a theory of nucleons. The justification for this is that we 
consider momentum transfers too small to resolve the structure of an individual nucleon. This 
matching is accomplished by considering the matrix element of the QCD operator between 
nucleon states. 

For the scalar operators qq, the result is well-known to be related to the so-called pion- 
nucleon sigma term [52, 53]. The nucleon mass arises from the trace anomaly: 

M n = (n|0» = Vm g (n\qq\n) + ^ {n\G%G a ^\n) , (4.6) 

q 4 ° s 

so one defines the fraction of nucleon mass arising from the various quark flavors and the 
gluons as: 

An) _ (n\m q qq\n) 



q=u,d,s 



Thus both the scalar quark operators qq and the scalar gluon operator TrG 2 match onto 
the scalar nucleon operator fin with coefficients that depend on how the various quark flavors 
contribute to the nucleon mass. A recent discussion of uncertainties in these quantities appears 
in [54]. 

The vector current case is easy to understand because we have a conserved current: 

qi^q ->■ Qnfij^n, (4.9) 

where Q n is the charge of the nucleon under the corresponding current. In particular, only 
valence quarks contribute at leading order in this case, not sea quarks or gluons. On the other 
hand, heavy-quark currents can lead to dipole moments even though the nucleon has no net 
charge under these currents; such dipole moments have been discussed in [55]. A subtlety 
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arises in computing them. From eq. 4.3 one might naively expect an m^ 4 suppression of the 
matrix element. However, the 3-gluon operator we matched to was renormalized at a scale 
uiq, and its matrix element in the nucleon has a quadratic divergence, so that the dipole 
moment is in fact suppressed only by irig 2 . 

The next familiar case is the axial-vector current qj^^q, which is related to the nucleon 
spin [40]. The matching is: 

<Z7^75<? -> A 9 n7 M 7 5 n, (4.10) 

where A q is the fraction of the spin of the nucleon carried by quark flavor q. These fractions 
are given by integrals of helicity-dependent parton distributions, as reviewed in [56], which 
are well-measured [57, 58, 59]. The case of the pseudoscalar current q^q is related by PCAC 
(or a Goldberger-Treiman relation) to that of the axial current [53]. For the axial current of a 
heavy quark, we match to the operator involving gluons, as in eq. 4.4. The matrix element of 
these gluonic operators has been argued to be related to a higher-twist matrix element that 
is known (and small) [60]. 

The tensor operator qa^ v q is of interest if the nucleus has a dipole interaction with some 
new gauge boson. This operator has appeared in the DM literature, often in the context 
of a general operator analysis (including, but not necessarily limited to, refs. [61, 62, 63]), 
but to the best of our knowledge the only correct discussion in the DM context of its matrix 
element in the nucleon is in ref. [64]. This operator is familiar in the spin physics context, 
and is related to the "transversity distribution" h\(x) of the nucleon [65, 66]. In particular, 
although the NR limit of the operators ^7^75^ and ipcr^ip agree, in the relativistic limit 
they measure different quantities, and the matching of quark-level operators to nucleon-level 
operators is different in the two cases. The analogue of the spin fraction in eq. 4.10 is the 
"tensor charge" 5 q : 

qo^vl^q -» Sqfia^^n. (4.11) 

Despite being twist-two parton distribution functions on an even theoretical footing with the 
well-studied cases, the transversity distributions are more difficult to measure and relatively 
poorly known. The tensor charges for light quarks in the proton have only recently been 
extracted from experimental data [67], giving, at Q 2 = 0.8 GeV 2 , 

5 U = 0.541°;°° (4.12) 
5 d = -0.23t°;°6- (4.13) 

4.3 Nuclear form factors for nucleonic currents 

Now that we have matched our effective theory onto operators involving nucleons, it remains 
to evaluate the correlators of nucleonic currents in nuclei. For instance, we may wish to 
calculate (N(p + q)\nn(q)\N(p)) . It is at this stage that we can take an NR limit and view 
this as simply a distribution of the number of nucleons or the spin of nucleons within a 
nucleus. These are familiar form factors [43, 45]; see refs. [28, 29] for more recent discussions 
of the status of our knowledge of the SI and SD form factors in various isotopes. 
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In plotting limits we have used the usual SI form factor F 2 (Er) and SD form factor 
S(Er)/S(0). In general, however, the story is more elaborate. Note that at this stage of the 
matching, multiple form factors can arise for a given process, e.g. the familiar decomposition 
for a coupling to a conserved current: 

(N(p + q)\n<fn(q)\N(p)) = N(p + q) (f^ 2 ) {2P + qY + {F^q 2 ) + F 2 {q 2 )) N{p ). 

\ IrriN rriN J 

(4.14) 

Here we have suppressed Lorentz indices on the spin- J nucleus wavefunctions N(p), which 
have the relativistic normalization (e.g. um{p) for spin-1/2 nuclei), and T,^ u is the spin-J 
Lorentz representation (e.g. gcr^ for fermions). It is only at this stage, when we have an 
effective coupling to the entire nucleus, that we can take NR limits of the spinor structure, 
finding that the 7* term is suppressed by the nucleus velocity relative to the 7 term measuring 
the total charge. In particular, note that the nucleon velocity v n ~ 0.1 does not appear, as it 
is not possible to take an NR limit at an earlier stage of the matching process. 

The second, magnetic dipole moment, term is suppressed by q/m^ ~ 10~ 3 relative to 
the leading term, assuming F2{q 2 ) is not much larger than F\{q 2 ). However, depending on 
which DM operator cpf^q is contracted with, it is possible that the two terms contribute in 
comparable ways. This possibility has recently played a role in several very recent attempts 
to understand whether data from DAMA, CoGeNT, and other experiments are compatible 
[10, 11, 12, 13]. If the coupling is to xo iliV \q h ', then in the NR limit the dominant part of the 
DM operator is ~ (s x x q), which is dotted into the velocity suppressed part of the quark 
current ~ ZvF\(q 2 ). On the other hand, the 7 term in the quark current contract with a 
velocity suppressed part of the DM operator. Thus F\{q 2 ) multiplies the operator s x -(vxq) in 
the NR limit. However, the magnetic moment term can be contracted without extra velocity 
suppression, corresponding to the operator F2(q 2 )q 2 s x ■ s/v in the NR limit. In this way, the 
F\{q 2 ) and F2(q 2 ) contributions can be comparably important in certain models. 

5. Conclusion and outlook 

In this paper we have constructed a simple NR effective theory to study signals from DM 
direct detection experiments. Different operators in the NR theory correspond to different 
types of interactions between DM and SM quarks. They lead to qualitatively different recoil 
spectra. Thus if DM is discovered in the near future DM direct detection, the recoil spectrum 
will constrain the NR effective theory and its possible field theory completions. Valuable 
information on the nature of DM-nucleus interactions can be obtained. 

There are still many details that need to be studied to optimize our understanding of 
possible DM signals. For instance, both the nuclear form factor and the DM velocity distri- 
bution could modify the recoil spectrum and smear the differences from DM dynamics. As 
far as we know, a complete list of the matching between quark or gluon operators and nucleon 
operators and the nuclear form factors is still absent in the literature. Another issue comes 
from the DM velocity distribution. In this paper we used the traditional simplest Maxwellian 
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distribution. But the DM velocity distribution in our galaxy is certainly more complicated, 
and non-Maxwellian distributions may result in more structures in the recoil spectrum. It 
will be interesting to develop methods to disentangle these nuclear and astrophysical effects 
from dynamics in analyzing data. (Some recent papers emphasizing the important role of the 
velocity distribution, among other uncertainties, were refs. [68, 69].) 

Finally, the interaction that triggers signals in DM direct detection may leave imprints 
in hadron colliders like the Tevatron and LHC. The effective theory differs for these two 
types of experiments due to the different kinematic regimes that they probe. Still, for certain 
interactions, there may exist strong correlations between the direct detection and collider 
signals. Then it will be highly desirable to extract and compare information on DM dynamics 
from both classes of experiments. 
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A. Direct detection recoil rate 

The basic quantity that direct detection (except for bubble chamber type experiment which 
is mainly senstive to the total rate) measures is the differential event rate per unit recoiling 
energy 



Nt is the number of scattering centers per unit detector mass, v is the DM speed relative to 
the target nucleus and p x is the local halo density of DM particle near the Earth, which we 
take to be p x « 0.3 GeV/cm 3 . Then vn x = vp x /m x is the incident WIMP flux, and f(v , ve) 
is the WIMP velocity distribution in the galactic halo with velocity in the range 500km/s < 
ve < 600km/s. In the numerical studies, we take the distribution to be Maxwellian with 
ve = 500 km/s, and average over modulation due to the solar system and Earth velocities 
[45, 70] during the dates the experiment ran, when reported, or over the whole year when not 
reported. The minimum velocity to scatter with recoiling energy Er is 




(A.l) 



Vmin — 



y/2m N E R 
2p N 



+ 



s/2m N E R 



8 



(A.2) 
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where is the nucleus mass; /xjv is the DM-nucleus reduced mass and for the inelastic DM 
scenario, 5, usually a mass splitting between DM components, characterizes the energy lost 
in the scattering to the dark sector. The total number of counts then follows from integrating 
over the energy bins, and multiplying by the effective exposure, e.g. in kg-days, for a given 
detector. Throughout the paper, we focus on elastic scattering and set 5 = 0. 

The form of the counting rate (A.l) isolates the physics of the WIMP-nucleus interaction 
in the differential cross section da /dER. This differential cross section is computed from the 
NR matrix element as (averaging and summing over initial- and final-state spins): 



da 

dER 



m N 



1 



2ttv 2 (2J + 1)(2s x + 1) 



NR\ 



(A.3) 



spins 



We have denoted the spin of the nucleus as J. As mentioned in section 2.2, for concreteness 
we will specialize to the case that x is a Dirac fermion, so s x = 1/2. 

Combining Eq. (2.1), Eq. (2.2), Eq. (2.3), we have the recoil rate in terms of the param- 
eters of NR effective theory as: 



dR 

dER 

da SI 
dER 

da SD 
dE R 



N T 



m 

\2 772 



L'E 



d 3 vvf(v,v E ) 



da 

dER 



A'F\E R )m N 



2irv 2 

J(J + l)m N S(E R 



h x + 



2irv 2 



5(0) 




h2\j2m^E R + 



y/2m N E R 
h' 2 \/2m N ER + 



y/2m N E R 



(A.4) 



Note that the terms in V e ff do not all interfere with one another, due to their differing spin 
structure. 



B. Power counting 

In this Appendix, we do some simple power counting exercises for the two limits of mediators 
to show that the leading momentum suppressed operators could still be relevant to direct 
detection. As discussed in Sec. 2.1, there are at most four scales relevant to DM scattering 
off nucleus: DM mass m x , nucleus mass mjv, mediator mass and some hidden sector 
scale A' unrelated to the first three. The NR expansion parameters fall into two classes: 
v, q/m x , qjvnjs! ~ 10 -3 for DM mass around the weak scale, and q/m^, q/A' which are unfixed. 
As argued in the context of momentum dependent DM, q/m^ could be as large as 0.1 as in 
the GeV dark sector models. Notice that as we argued in Sec. 4, there could not be large 
expansion parameters such as the nucleon velocity v n ~ 0.1 present. 

For elastic SI scattering, the smallest DM-nucleon cross section probed by the present 
direct detection experiments is around a p ~ (10~ 43 — 10~ 44 ) cm 2 . Future upgraded experi- 
ments such as XenonlT will probe cross sections down to 10 cm 2 , close to 10 8 cm 2 , the 
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irreducible neutrino-induced nucleus recoil background. For elastic SD scattering, the current 
bound on cross section is around the same order of magnitude of that of the weak interac- 
tions, a p ~ (1CT 37 — 10 -38 ) cm 2 . (Despite the orders of magnitude separating these limits, it 
is not at all unnatural to expect both SI and SD signals; for instance, this is a very robust 
expectation if the DM is the MSSM LSP [71].) The bounds on the averaged DM-nucleus 
cross section are 




where we defined the averaged DM-nucleus cross section to get rid of the multiplicity factor 
of atomic number A or total charge Z. 

Now we do a crude power counting exercise to estimate to what order the operators could 
still contribute to direct detection: 

• Heavy mediator: ~ 0(100 GeV). 

iMso ~ h'l, (B.3) 

where the dimensionless coefficient h{h!) = g x f nuc \ represents products of DM and 
nucleon coupling to the mediator. From Eq. (B.l), we see that to have experimentally 
accessible cross sections, 

hi ~ 10" 4 , h 2 ~ 0.1, h[ ~ 0.1. (B.4) 



Light mediator: <C q 2 . 



iMsi ~ \h + ■ ■ ■ + kO i^— 
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where A collectively denotes DM mass or nucleus mass around 100 GeV. Here the 
couplings l\ ~ 10~ 10 , ~ 10 3j /i, 1' ~ 10~ 7 , ~ 10 31 l' Q are relevant for direct 



In both cases we truncate the expansion when the coefficient has to be of order one to be 
relevant for direct detection. More accurate bounds on the coefficients have been presented 
in Sec. 3.1. 

C. Matching to effective field theory operators 

In this appendix, we will discuss the general rules to write down NR effective operators and 
the matching between the effective field theory operators for DM with different spins and the 
NR effective theory. 

C.l Fermionic DM 

C.l.l Complete set of NR operators 

In the center of mass frame, the amplitude of fermionic DM scattering off a nucleus can 
be expressed in terms of scalar invariants formed out of four independent 3-vectors: the 
transferred momentum q, the DM incoming velocity v relative to the target and the DM and 
nucleus spins s x and sjy- In fact, it is somewhat more natural to parametrize the operators 
not in terms of v but instead in terms of the sum of momenta before and after scattering, 
pt 1 = (p -\- p'^ = (2p + q)^, which is related to the velocity by P = 2[inv + q. With 
two spins and two momenta, it is found in [72] that 16 independent rotationally invariant 
operators can be constructed and they include all possible spin configurations. Any other 
scalar operator involving at least one of the two spins can be expressed as a linear combination 
of the 16 operators with spin-independent coefficients that may depend on the momenta 
through \q\ 2 , \ v\ 2 . Below is the list of the complete set in the momentum space without 
mediator propagators, with superscripts denoting P and C parity. In particular, notice that 
there are six CP- violating operators: O2, Oq, O9, O13, O14, and O15. 

• SI operators 



detection. 



O 



(++) 



= 1 



1 




(-+) 



= is x ■ q 

= s x -P 

= is x ■ (P x q) 



(C.l) 
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SD operators 



-.(++) 


= s x -s N 
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= is N -q 







0<-> 


= s N -P 







{—) 


= i(s x x s N ) 


• q 
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= (s x x s N ) 
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= isn ■ (P x 
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11 = {s N -q)(s x -q) 
,{ + 2 +) = (s N -P)(s x -P) 

,(+-) 

13 

,(+-) 
14 

,(-+) 
15 

)(— ) 
16 



(sjy ■ q)(s x ■ P) + (sat • P)(s x 
(s N ■ q)(s x ■ P) - (s N ■ P)(s x 



s N -(Px c 
s N • (P x q) 



1) + 
P) + 



(P x 



(sn ■ q) 

Pxq)^ (s N ■ P) 
(C.2) 



Among the 16 operators, there are four nucleus SI operators {0\ — O4) with only O4 
suppressed by two powers of the expansion parameter \v\, \q\. In coordinate space, 0\ — O3 
lead to 3 types of contact interaction for heavy mediator and 3 long-range interactions for 
light mediator. The 4 cases arising from the two operators 0±, O2 (with at most a single |<f| 
suppression) are included in our parametrization in Eq. (2.2). 

The other 13 operators O5 — 0\% are all SD. Up to leading order in v, q, one has O5 — 
Og . The amplitudes of Os , O9 are almost the same as those of Oq , O7 except for different 
dependence on the DM spin. Thus, after averaging over initial spins and summing over final 
spins, 08(0g) give differential cross sections identical to those from 06(07) up to a numerical 
factor. Thus we do not include them in Eq. (2.3). 

The physical interpretations of the operators are more transparent in the potential form. 
For instance, if the mediator is light, 02, Oq are the dipole-monopole potentials while On is 
the dipole-dipole potential. 



C.1.2 Field theory operators 

Now we turn to field theory operators and write down the most general four fermion interac- 
tion mediated by some scalar or gauge boson in the form: 

G(q 2 , m^)J x J q = G(q 2 , m <t> )xT x xqT q q 1 (C.3) 

where J x and J q are the DM fermion and quark bilinears (we hope the reader will not 
be confused by our use of q both for quarks and the momentum transfer). The propagator 
G(q 2 ,m x ) is approximated as l/m^L in the heavy mediator case and 1/q 2 in the light mediator 
case. The DM fermion could be a single Majorana, Dirac or two Majoranas. is a 4 by 

4 matrix, which we choose to be 1, 7 s , 7^, 7^7 5 , , a^ u j 5 for non-derivative interactions. 
Correspondingly, we have the following DM and quark non-derivative interactions 

(pseudo)scalar x (pseudo) scalar : x{l^)xq{lb)q 
(pseudo) vector x (pseudo)vector : Xl IJ '{l b )xqi^i.l 5 )q 
(pseudo)tensor x (pseudo)tensor : x^ 1 '(j 5 )xq~° 'tiu(l 5 )q- 
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Effective operator 


leading NR operators 
in momentum space 


leading NR operators 
in position space 


SI 


xxqq 
ixi 5 xqq 

xi 5 i^xqifiq 

ixa^^xqcr^q 


1 

is x -q 
1 

is x -q 


5 3 (r) 
-s x • W 3 (r) 
5 3 (r) 

^x^+^s x -V)5 3 (r) 
-s x • W 3 (r) 


l 

r 

s x -r 

r a 
1 
r 

s x -v iS x -r 
r 2^tjvr 3 
s x -r 
r 3 


SD 


ixxqisq 
xi b xqi b q 

xi^xqi^q 
xr'i^xqi^q 

Xa^xqv^q 


is N ■ q 

(s x • q)(sN ■ q) 

i(s x x sn) • q 

s N -v 
i(s x x sn) • q 

s x ■ s N 

s x -s N 


-s N ■ W 3 (r) 
(s x -V)(s N -V)5 3 (r) 
-(s x x s N ) • W 3 (r) 

-(s x x s N ) ■ W 3 (r) 
s x ■ s N 5 3 (r) 
s x ■ s N 5 3 (r) 


3(s x -r)(s]v-r0 

r-(s x Xs N ) 
47rr 3 

f-(s x Xs N ) 
47rr 3 

r - 
r 



Table 2: Effective operators of fermion DM and their corresponding NR operators 



Notice that not all of them are independent. For instance, as = |e^ pc 'a pa , the pseudo- 

tensor x pseudo-tensor interaction is equivalent to the tensor x tensor interaction and only 
one pseudo-tensor x tensor interaction is independent. If the DM is a single Majorana 
fermion, operators such as xi^X an d X®^ '{l 5 )x do not exist. But in the NR limit, one 
cannot tell the difference between single Majorana and Dirac fermions as they lead to the 
same set of NR operators. Below we summarize in Table. 2 the correspondence between field 
theory operators and NR operators. We also list both contact and long-range potentials for 
every operator. 

A few comments are in order: 

• We exhausted all possible four fermion non-derivative interactions, but the NR leading 
operators in the momentum space only cover 8 of the complete set of the 16 operators 
listed in Eq. (C.l, C.2). In order to get the other operators, one needs to combine and 
manipulate the effective operators listed above so that the NR leading terms cancel out. 
For instance, to get the SI operator O4 = is x ■ (P x q), tuning is necessary to cancel 
a leading operator. (Note that, e.g., an expression that Fierzes to xxqq ~ Xl^XTI^q 
does not lead to an automatic cancellation, because the first term depends on A times 
nuclear matrix elements whereas the second depends on Z times different nuclear matrix 
elements.). 

• As noted above, the velocity appears through the combination P = 2fiNV + q, with 
UN the DM-nucleus reduced mass. For instance, consider a single Majorana fermion 
coupling to the Z boson. After integrating out the Z, the leading unsuppressed op- 
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erator is SD, X7 7^X(Z7 7^?> giving a potential s x ■ s/vc) 3 (r). The leading SI operator 
is X7 5 7 aj X97m9i giving V ~ (s x ■ v + j^^x ' V)<5 3 (r) ~ • v5 3 (f) where s 1 - is the 
spin transverse to the momentum transfer. This is known as DM "anapole" moment 
scattering off the nuclear electric current. The SI scattering cross section is suppressed 
by (q/mz) 2 ~ 1CT 6 . But notice that the current SD bound from direct detection is 
about 10 -38 cm 2 , six orders of magnitude above the SI bound 10 -44 cm 2 . Thus with 
all couplings of order one, both operators are relevant for direct detection. This is an 
example where one cannot naively discard a momentum-suppressed SI operator if the 
leading operator is SD. 

• Derivative effective field operators could also give NR operators in Eq. (2.2, 2.3). These 
operators arise naturally from dark moment models. For instance, the DM electric 
dipole moment scattering off the nucleus charge is encoded in the matrix element 

In the NR limit, it yields the (DM) dipole coupling s x -r/r 3 . The DM magnetic moment 
scatters off both the nucleus charge and spin, leading to NR operators s x ■ s/v and 
s x ■ (v x q) as discussed in the end of Sec. 4.3. 

C.1.3 Examples of matching 
Pseudoscalar Exchange 

We consider a light pseudoscalar coupled to the DM and visible sector as 

iac/)X7 5 X + P&qq + h.c, (C.5) 

where we assume that j3 is flavor diagonal. The coupling to quarks could be translated to 
that to the proton as 

q=u,d,s H \ q=u,d,s I q=c,b,t 

Then the coupling to the nucleus is /jv = Zf p + (A — Z)f n with form factor set aside for the 
moment. 

In the NR limit, the matrix element is 

M SI = iaf N 4m N q- s x . (C.7) 

After taking into account different normalizations, the NR matrix element is related to the 
relativistic matrix element as 



The corresponding NR effective potential is the Fourier transform of the NR matrix element 
to the leading order of Born approximation 



V SI = i 2 Sx 



4 7rr 3 

a fZ Z- A 



where l 2 = — ( -f p + -^—f n ) . (C.9) 



The differential recoil rate is then 



da SI _ A 2 F 2 {E R )m N 
dE R ~ 8^ 



y/2m N E R 



(CIO) 



DM anapole moment scattering 

Now we will move to a more elaborate example of the NR effective theory, which goes beyond 
the small set of operators that we considered in the main text. The model we consider has 
a Majorana fermion as DM. It interacts with the nucleus through the Z boson. This leads 
to the anapole moment scattering which has been considered recently, for example, in [12]. 
After integrating out the Z, the effective theory operator for the DM-quark interaction is 

ax7 M 75X<rrM<? + Pxi^ibxqi^im, (di) 

where the coefficients a, (5 are 

g^G 1 / (1 4 . 2 

an = — — ---sin 9 W 

mz \2 3 

9j ^2Gf ( 12 . 2 



mz V 2 3 

9x" 



,^2G]l 2 



2m z 

? d = ( c - 12 ) 

2m z 

where the subscripts of a, f3 denote the couplings of up or down type quarks. g x , the DM- 
Z coupling constant, is a free parameter. (As mentioned in Sec. 2.3, it can arise from a 
higher- dimension operator such as xi^X^ D^h, and so is not directly related to the SM 
gauge couplings g and g' .) Then we use the recipe outlined in Sec. 4 to map the DM-quark 
operator to DM-nucleus operator. At the level of nucleons, we have couplings: 

X7^75X (Jfi + AO = X7 M 75X (fpPlfiP + funj^n + a p pj^ 5 p + a„n7 / ,7 5 n) , (C.13) 

Where we have denoted the vector and axial currents on the right-hand side by J and A and 
defined: 

f P = 2a> u + a d , 
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fn = oc u + 2a d , 
a p = 



q = u,i 



g = u,d,s 



(p) 

with Ag the fraction of spin of the proton carried by quark flavor q, as reviewed in Sec. 4.2. 
Now, we evaluate matrix elements in the nucleus: 

X1^X{N\J,\N) = ZU+ ^~ A Z)fn Xl^xN [AF(q)(p + p') ii + 2i^F( q ))N 

X7 M 75X (N\A,\N) = (a p ^ + aj-^j xi^xN{2S,)N. (C.15) 

Here we have defined AF{q) and F(q) to be the charge and magnetic dipole moment form 
factors for the baryon number current (we have assumed that protons and neutrons are 
distributed uniformly throughout the nucleus). denotes the nucleus spin. As in Sec. 4.3, 
N and N denote appropriate spin-J wave functions and a spin- J Lorentz representation. 
The two terms arising from J do not mix with each other or with the term arising from A, 
so will discuss the three separately. 

For convenience, we will define constants: 

Zf + (A - Z)f n 
C J = — 



A 

(S p ) (S n ) . . 

ca = a p —j- + a n —j- . (C.16) 

Now we will calculate the NR matrix elements. First, from J we extract the piece proportional 
to F(q). A straightforward calculation shows that it gives a spin-independent matrix element. 
Here we factor out AF(q), according to our convention for SI NR matrix elements as set out 
in Sec. 2.1. (Recall that we have assumed x to be Majorana; certain factors of 2 would differ 
if it were Dirac.) 

i SI 



M = 16m x m,NCj I v + 



2hn, 

= 16m x mNCj(s x — q ■ s x q) ■ v = 16m x mNCjs^ ■ v. (C17) 

Next, we consider the spin-dependent piece arising from the magnetic dipole moment term in 
the matrix element of J . Because it is spin-dependent, we follow our convention of factoring 
out \j even though the form factor in this case, F(q), is different. This leads to some 
awkwardness in the definition of our coupling: 



M S j D = McjF(q)J j^m x s x ■ (qx s N ) . (C.18) 
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Notice that this matrix element corresponds to the operator 0% in the general classification 
of App. C.l.l. Finally, we have the spin-dependent matrix element arising from the axial 
current A: 



M S j D 



(C.19) 



Then, dividing by Am x m^ to convert to the nonrelativistic normalization, we obtain the 
corresponding NR effective potentials: 



V 



si 



h i ] v-s x + ih 2 V ■ s x ) 5 3 (r), 



V 



SD 



h\ s N ■ s x 5\r) + h 8 (s x x s N ) • V5\r), 



(C.20) 



where 



k?=2 mh2 = 4^^ A A - Z ^ 

h ' l = 8 (^A + a J s - 



j " j 

Zf p + (A — Z)f n 



S(q)' 



(C.21) 



The differential recoil rate is then 



da 



si 



dE R 

da SD 
dE R 



A'F\E R )m N 

2lTV 2 



\ h i )y + ]jh 2 \/2m N E R 



J{J+l)m N S{E R ) (I 



2ttv 2 



S(0) 



4 1 11 12 



h 8 \/2m N E R 



(C.22) 



Notice that the SI rate is slightly different from Eq. (A. 4) as it includes the interference term 
between the velocity-suppressed and momentum-suppressed operators. Furthermore, the SD 
rate includes a term very similar to the familiar h' 2 contribution, but arising instead from the 
operator Og = i(s x x sjv) ■ q- Thus, this example goes beyond the formalism of App. A. We 
have avoided such subtleties in most of the text because the limited set of coefficients hi, 
h'i, li, l\ suffice to produce most of the interesting variations in recoil spectrum shape. Our 
new hi term gives a shape similar to that arising from hi, while h% gives a shape similar 
to that arising from h' 2 . The presence of these new coefficients in an example as natural 
as Z exchange, however, shows that the inverse problem of matching a field theory to a 
measured shape is more subtle, and necessarily involves a number of degeneracies due to the 
higher-dimensional space of possible potentials that yield similar shapes. 

C.2 Scalar DM 

The analysis for the scalar DM is greatly simplified due to the trivial fact that scalar has no 
spin. In the momentum space, only three NR operators could be constructed from q, v, sjv to 
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the order of \q\, \ v\ : 

Osi = 1 

Osd = s N ■ q + s N ■ P 

As we only consider static potential, only the first two operators remain which, in the position 
space, are translated to 



Vsi = 5\r), 



47rr 



V SD = s N ■ V<5 3 (r), 



isn ■ r 
4irr 3 



(C.23) 



In the field theory, one could write down scalar operators that couple to the (axial) scalar 
and (axial) vector quark current: 



They are related to the NR operators in Table 3. 



(C.24) 





Effective operator 


leading NR operator 
in momentum space 


leading NR operators 
in position space 


SI 


4> 2 qq, ftd^qj^q 


1 


6 3 (r) 


1 


SD 


<t> 2 qj 5 q 
tfd^qj^q 


sn ■ q 
s N -v 


sn ■ V5 3 (r) 
(sN-v + ^s N -V)5 3 (f) 


isfff 

47C7-3 
s N -v is N -f 



Table 3: Effective operators of scalar DM and their corresponding NR operators. 



C.3 Vector DM 

The spin-1 DM NR operator analysis is analogous to the fermionic DM. One modification is 
the representation of the spin operator. In the NR limit, vector DM spin operator is 

sb = e ' • R ■ e 

or s| = e ijk e*ej (C.25) 

where e % is the 3 component polarization vector and R is 3-dimensional spin-1 rotation repre- 
senation with (R k )ij = e l ^ k . In addition, besides the vector representation, high-dimensional 
tensor representation of the rotation group could be relevant (one do not need to consider 
this for the spin 1/2 case as the Pauli matrices form a complete basis for the 2 by 2 matrice). 
For instance, as we will show below, the symmetric representation with two indices could be 
relevant 

S ij = J ■ R {i R j) ■ e = e+M . (C.26) 
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One prototype example of vector DM in the field theory is the the first KK mode of the 
photon in the universal extra dimension model [73]. Again, one can write down operators 
involving the heavy vector DM BP that couple to the scalar and vector quark current. First 
if B^ is self-conjugate, 

BPB^)q, 

BPdpB„qY(7 5 )q ,epv P aB»dPB°q 7 P( 7 5 )q. (C.27) 

Additional operators exist if the gauge field is complex: 

B^d v B^tf)q, (C.28) 

These operators are not gauge invariant and are generated from UV theory after integrating 
out heavy degrees of freedom. Thus we only consider contact interactions. In the NR limit, 
they lead to operators in table 4. 





Effective operator 


leading NR operator 
in momentum space 


leading NR operator 
in position space 


SI 


BPB^qq,B^3 v BpqYq 
e^B^dPB^q^q 

BPdpB u qYq 


1 

s^-v 
iP-S-q 


5 3 {r) 

(sW+^ B -V)* 3 (r) 
-(^+5^V)-5-W 3 (f) 


SD 


BPB^qfq 
B^d u BpqYl b q 

e llvpa B v dPB°q 1 P 1 5 q 
er,p a B v dPB a vfq 


iq- s N 
s N -v 
isN • S • q 
s B ■ s N 
i(s B x s N ) -q 


-s N ■ V5 3 (r) 
(sN-v + ^s N -V)5 3 (f) 

-8 N -S- V5 3 (f) 

s B ■ s N 5 3 (r) 
-(s B x s N ) ■ <5 3 (r) 



Table 4: Effective operators of vector DM and their corresponding NR operators. 



We will not explicitly consider higher-spin DM, although the story generalizes straight- 
forwardly. A discussion of the spin-3/2 case can be found in [74]. 
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